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ABSTRACT 

We report on a BeppoSAX observation of the transient X-ray pulsar 4U 0115+63 close to periastron. 
This led to the discovery of a dramatic luminosity variation from ~ 2 x 10'^^ ergs~^ to ^ 5 x 10'^^ ergs~^ 
(factor ^ 250) in less than 15 hr. The variation was accompanied by only minor (if any) changes in 
the emitted spectrum and pulse fraction. On the contrary an observation near apastron detected the 
source in a nearly constant state at a level of '-^ 2 x lO'^^ ergs~^. Direct accretion onto the neutron star 
surface encounters major difhculties in explaining the source variability properties. When the different 
regimes expected for a rotating magnetic neutron star subject to a variable inflow of matter from its 
companion are taken into consideration, the results of BeppoSAX observations of 4U 0115+63 can be 
explained naturally. In particular close to apastron, the regime of centrifugal inhibition of accretion 
applies, whereas the dramatic source flux variability observed close to periastron is readily interpreted as 
the transition regime between direct neutron star accretion and the propeller regime. In this centrifugal 
transition regime small variations of the mass inflow rate give rise to very large luminosity variations. 
We present a simple model for this transition, which we successfully apply to the X-ray flux and pulse 
fraction variations measured by BeppoSAX. 

Subject headings: stars: individual: 4U 0115+63 — stars: neutron stars — X-ray: stars 



HXRT emerges out of quiescence or fades away at the end 
of an outburst therefore holds a great potential as a diag- 
nostic of the different regimes experienced by a rotating 
magnetic neutron star. 

In this paper we report on BeppoSAX observations of 
the transient X-ray pulsar 4U 0115+63, which reveal the 
quiescent state of the source and, more crucially, provide 
the first convincing evidence for the centrifugal transition 
regime of any X-ray transient. 

2. 4U 0115+63 

Like most HXRTs, 4U 0115+63 hosts a magnetic neu- 
tron star orbiting a Be star in a moderately eccentric orbit 
(e ~ 0.3). The spin and orbital periods are Pq = 3.62 s 
and Porb — 24.3 d (Cominsky et al. 1978; Rappaport et al. 
1979). The system has been extensively studied during its 
frequent outbursts {^ 20 have been detected so far) , which 
occasionally reach a peak luminosity of ~ lO'^^ ergs""'^ 
and usually last about a month (Bildsten et al. 1997; 
Campana 1996). The neutron star magnetic fleld is mea- 
sured with good accuracy, Pq = 1-3 x 10^^ [(1 + z)/1.2] 
G, through the cyclotron features detected in its X-ray 
spectrum (Santangelo et al. 1999 and references therein). 
Here (1 + z) = l^GM/{c^ R) ~ 1.2 (we scale the neutron 
star mass and radius as M = 1.4 Mi. 4 Mq and R = IORq 
km, respectively) is the gravitational redshift at the neu- 
tron star surface. The source distance was determined to 
be (8 ± 1) dg kpc, such that the X-ray flux to luminos- 
ity conversion is correspondingly accurate (Negueruela & 

^ Osservatorio astronomico di Brera, Via Bianchi 46, 1-23807 Merate (LC), Italy 

^Dipartimento di Fisica G. Occhialini, Universita di Milano Bicocca, Piazza della Scienza 3, 1-20126 Milano, Italy 

^Osservatorio astronomico di Monteporzio Catone, Via Frascati 33, 1-00040 Monteporzio Catone (Roma), Italy 

4 Affiliated to I.C.R.A. 

^TeSRE, C.N.R., Via Gobetti 101, 1-40129 Bologna, Italy 

1 



1. INTRODUCTION 

Accreting collapsed stars in transient X-ray binaries are 
subject to very large variations of mass inflow rate and 
provide a laboratory to test the physics of accretion over 
a range of different regimes that is unaccessible to persis- 
tent sources (see e.g. Parmar et al. 1989; Campana et 
al. 1998). The highly magnetic (P -- lO^^ q)^ spinning 
neutron stars that are hosted in Hard X-ray Transients 
(HXRTs) have long been suspected to be in the so-called 
propeller regime when quiescent. In this regime accre- 
tion onto the neutron star surface is inhibited by the cen- 
trifugal action of the rotating magnetosphere (lUarionov 
& Sunyaev 1975). On the other hand there is little doubt 
that accretion onto the neutron star surface takes place 
in these systems while in outburst, as they display very 
similar properties to those of persistent X-ray pulsars. 

Relatively little is known on the low luminosity states 
of HXRTs. Observations of their quiescent state are still 
sparse (Campana 1996; Negueruela et al. 2000) and there 
are only hints that the transition from the lowest outburst 
luminosities to quiescence occurs in a sudden fashion (see 
e.g. the case of V 0332+53, Stella et al. 1986). Accord- 
ing to models of accretion onto rotating magnetic neutron 
stars (e.g. Illarionov & Sunyaev 1975), the transition from 
the accretion to the propeller regime (or vice versa) should 
take place over a very limited range of mass inflow rates 
and give rise to luminosity variation in the order ^ 1000 
for neutron star spinning at a few second period. The 
observation of a paroxysmal luminosity variation when a 



Table 1 
Exponential cut-off power law fits to spectra from different intervals from the August 3-4, 1999 observation. 



Interval 


Time 


Nh 


Power law 


-E-cutoff 


-Efold 


0.1-200 keV Flux!i 


-^rcd 




(hr) 


(1022 cm-2) 


r 


(keV) 


(keV) 


(erg s"-*- cm~2:) 


(d.o.f.) 


1+ 


7.5-14.9 


2 30+ -^■"** 


1 t;«+U-4U 
J^-O»_0.40 


12.0 (fixed) 


8.2 (fixed) 


1.1 X 10"" 


1.52 (24) 


2 


15.4-17.2 


l-83_();57 


0.82tO-i8 


14.3l2;| 


6.3l5;2 


1.5 X 10-10 


1.13 (28) 


3 


17.5-19.4 


1-80-0.64 


i.UU_Q 20 


11.71^2^0 


I3.2+1O35 


1.5 X 10-10 


1.07 (33) 


4 


20.5-21.5 


-I 7^+0.40 


o.9ot°:l3 


9.71?- 


8.5+2-3 


2.8 X 10-10 


1.04 (43) 


5 


22.0-22.5 


1 41+0-35 
^■^^-0.30 


0.89l°:l? 


10.31^:° 


s-otli 


5.0 X 10-10 


0.97 (40) 


6t 


22.5-22.7 


1.74 (fixed) 


1.101013 


12.51^3 


6.0l^-« 


5.7 X 10-10 


0.78 (41) 


7t 


22.7-23.3 


1.74 (fixed) 


0-95lS:w 


13.2+f* 


5.5+Jl 


7.3 X 10-10 


1.48 (36) 


8 


23.3-23.9 


1 OO+0.46 
^■°3_o.36 


0.94iril 


9.4it:^ 


ii.2li:i 


1.1 X 10-** 


1.02 (33) 


9 


23.9-24.2 


J^-^»_0.39 


o.s5+°:\l 


11.712-] 


7.Stli 


8.9 X 10-10 


0.98 (44) 


lot 


24.2-24.4 


1.74 (fixed) 


0.9310:11 


9-ltl^ 


9-7tl:t 


9.0 X 10-10 


0.88 (44) 


11 


24.5-25.7 


o 4^ + 1.00 


1.3010-21 


15.0l^,.o 


9 4+4.2 
y-^-3.1 


7.3 X 10-10 


1.13 (36) 


Mean value 




1.74 ±0.18 


0.95 ±0.05 


12.0 ±1.3 


8.2 ±1.1 






xtd (cons) 




1.7 


2.1 


0.8 


0.6 







Time is measured in hours from the start of MJD 51393. 

§ Fluxes are unabsorbed. Errors are at 90% confidence level (c.l.) for one parameter of interest (i.e. Ax^ = 2.71). 

+ Due to poor statistics, we fixed the cut-off and folding energies to the mean value of the entire observation. In the case of a simple power 

law fit, a column density of (2.2lJ-g) X 1022 cm-2 ^^^^ ^ photon index of F = 1. 541^-33 (x^ed ~ 1-44) were obtained. 

t Only MEGS and PDS data available. The column density has been fixed to the mean value since without the LEGS data it is difficult to 

constrain its value. 



Okazaki 2001). 

Before the present study, 4U 0115+63 was not detected 
in its quiescent state to a limit of ^ 6 x 10^'^ ergs~i (0.4-6 
keV; Campana 1996). At the end of the 1990 outburst a 
quick decrease of the luminosity was observed below a level 
of '^ 5 X 10^^ ergs"^, possibly indicating the transition to 
the onset of the centrifugal barrier (Tamura et al. 1992). 

3. BEPPOSAX OBSERVATIONS 

A 35 ks observation of 4U 0115+63 was carried out with 
the Itahan/Dutch satellite BeppoSAX (Boella et al. 1998) 
on 1999 August 3 06:37:04 UT - August 4 01:41:32 UT, 
covering an orbital phase range (0.94-0.97) close to perias- 
tron. The source count rate varied by a very large amount 
(from }^ 0.02 to ~ 5 c s~^) showing a steep increase by a 
factor of ^ 250 in ~ 15 hr (see Fig. ||; note that the Y- 
axis is logarithmic). Superposed to this trend, variations 
of up to a factor of ~ 20 in less than one hour were clearly 
present. 

Inspection of the RossiXTE ASM light curve of 4U 
0115+63 shows that the dramatic flux increase revealed by 
the BeppoSAX observation was not followed by a source 
outburst over the following days; the source remained in- 
stead below a daily average of ~ 2 x 10^^ ergs~^. 

In consideration of the extreme variability, we divided 
the entire observation in different time intervals (for spec- 
tral analysis) as well as intensity intervals (for temporal 
analysis^), in order to characterize the source properties 
at different luminosity levels. Energy spectra from the 
BeppoSAX LEGS (0.1-4 keV), MEGS (1.8-10 keV) and 
PDS (15-200 keV) were accumulated in each of these in- 
tervals. In all cases the source spectrum from the three 
instruments was well fit by a model consisting of an ab- 
sorbed power law with high-energy cut-off plus absorption, 
with photon index of F '-^ 1.0, cut-off energy £^cut ~ 12 
keV, e— folding energy Efdd ~ 8 keV. The column density 

^Pulsations were detected also at the lowest intensity interval, even 
this is why time intervals were used instead in the spectral analysis. 



amounts to Nh ^ 1-7 x 10^^ cm-^ and it is consistent 
with the galactic value. This model fits also nicely the 
spectrum of the entire observation and is compatible also 
with the source spectrum measured during outburst (see 
e.g. Nagase 1989). The best fit parameters are reported in 
Table 1 for each interval (see also Gastaldello et al. 2001, 
in preparation). We note that all spectral parameters are 
nearly constant throughout the observation (see Table 1). 
Variations are observed only in the column density (4% 
probability of getting a higher x^ by chance from a con- 
stant distribution) and in the power law index (probabil- 
ity of 1%). The main contribution to the variations of the 
power law index comes from the first interval. If we ex- 
clude this value, we obtain a probability of 10%. Marginal 
evidence was found for the second cyclotron harmonic at 
a centroid energy of 22.5 ± 2.5 keV in the spectrum from 
the entire observation (95% significance; Santangelo et al. 
1999 measured 24.16 ± 0.07 keV during outburst). Note 
that the first harmonic falls just in the gap between the 
MEGS and the PDS spectra. 

In this letter we estimate the source bolomctric lumi- 
nosities by using the 0.1-200 keV unabsorbed luminosities 
as derived from the best fit parameters of the spectrum 
from the entire observation and the MEGS count rates in 
the 1.8-10 keV band. The corresponding conversion factor 
was lcs~i = 8.5 x 10"^^ ergs~i. The luminosity inferred 
in this way varied between 2 x 10^'' and 5 x lO'^^ dg ergs~i 
(see Fig. g). 

Pulsations at the ~ 3.62 s neutron star spin period were 
detected in all intensity intervals. The pulsed fraction 
(semi-amplitude of modulation divided by the mean source 
count rate) in the 1.8-10 keV energy band was determined 
for each interval after folding the data at the best period. 
Values ranged between ~ 29% and ~ 54% with the lowest 
value corresponding to the lowest intensity interval (see 

Fig. D- 

though the source spectrum could not be meaningfully characterized; 
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Fig. 1. — Light curve of the BeppoSAX MECS observations of 4U 0115+63 taken in quiescence (left), during the transition (middle) and 
in outburst (right). Background subtracted light curves were converted into luminosities using the conversion factor derived from spectral 
fits (see text). Time bins are 10,000 s, 500 s and 5,000 s for the three observations, respectively. The time axis values (in hr) correspond to 
the start of the three observations discussed in the text (MJD 51769, 51393 and 51263, respectively). Solid lines represent the luminosity 
corresponding to the onset of the centrifugal barrier [Liim(i?)], and the luminosity at which it closes completely [Liim(^m)] for ^ = 1. Dashed 
lines give the same luminosities for { = 0.5. 



A further 86 ks BeppoSAX observation was carried out 
on 2000 August 13 21:33:06 UT - August 16 00:22:36 UT 
around apastron (orbital phase of 0.42-0.51). The source 
was very weak and remained at a virtuaUy constant level 
of ~ 2 X 10~^ c s~^ in the MECS; it was not detected 
in the LEGS (< 2 x 10"^ c s~^). The spectrum was 
softer than in the August 1999 observation. Fitting a 
power-law model to the MECS data and fixing the col- 
umn density to the best fit value of the periastron obser- 
vation yielded a photon index of T = 2.6^0 g (at 68% c.l., 
Xred ~ '-'•2 '^'-"^ 2 d.o.f.). The uncertainties in the spectral 
parameters translate also into a fairly large uncertainty in 
the inferred 0.1-200 keV unabsorbed luminosity, which is 
(0.6 — 3) X lO'^^ erg s~^. A pure blackbody model gave also 
a reasonable fit for a temperature oi kT = 0.7lg'2 keV 

(68%c.l.,x?ed = 0.5). 

The ^ 3.62 s pulsations were not detected; a 3 a upper 
limit of ~ 30% on the pulsed fraction was derived (see also 
Campana et al. 2001 in preparation). 

In order to facilitate the comparison with the source 
properties while in outburst we also analyzed the data 
from the 48 ks BeppoSAX observation that took place on 
1999 March 26, during the outburst decay. This observa- 
tion covered an orbital phase interval of 0.61-0.65. The 
source light curve is also shown in Fig. n for comparison; 
the average 0.1-200 keV luminosity was ~ 8 x 10'^^ erg s~^. 
The source spectrum could be described by a power law 
with a high energy cut-off with F ~ 0.8 ± 0.1, Ecut ~ 
9.4±0.5keV, £;foid ~ 16±2 keV and iV// ~ (1.5±0.1)xl022 



cm ^. To obtain an acceptable fit also an iron line and a 
cyclotron line feature had to be included. 



4. EVIDENCE FOR THE TRANSITION FROM THE 
PROPELLER TO THE ACCRETION REGIME 

The factor of ^ 250 luminosity variation of 4U 0115-1-63 
during the August 1999 close to periastron observation, 
when the source luminosity ranged between 2 x 10'^^ and 
5 X 10'^^ ergs~^ in ~ 15 hr is the most extreme ever seen 
in a HXRT. Previously the steepest flux increases detected 
from this source (and other HXRTs) were those associated 
to the rise to an outburst peak, involving variations of up 
to a factor of ^^ 3 — 4 on a comparable time scale. On 
the contrary during the August 2000 observation close to 
apastron the source luminosity was low and nearly con- 
stant around a level of 2 x 10^^ ergs~^. 

One possibility to explain the extreme variability in the 
August 1999 observation is that the source was subject to 
factor of ^ 250 variations in the mass inflow rate, which 
in turn gave rise to a comparable variation in accretion 
luminosity (under the hypothesis that accretion onto the 
neutron star surface took place unimpeded also for lumi- 
nosities as low as ^ 10'^'* ergs~^). This possibility faces 
serious problems. Firstly, models of Be star disks/winds 
predict a neutron star mass capture rate variation of a 
factor of ^ 5 at the most for a binary system such as 
4U 0115+63 over an orbital phase interval of 0.94-0.97 
(Raguzova & Lipunov 1998). Moreover, within the direct 



38 






o 



36 



34 



Llim(R) 




/ 




• 


• •• 




•\ 


.mT 












J^ 






Llim(rm) " 



10 



15 20 

Time (hr) 



25 






3 



0.8 


- 


- 


0.6 


1 


- 






^-+- 


^ ♦"-♦- 


0.4 


/ 
1 






-♦- 


0.2 


NK 1 












- 



33 



34 35 36 

Log Luminosity (erg s ) 



Fig. 2. — Left panel: Light curve of the 0.1—200 keV unabsorbed luminosity of 4U 0115+63 versus time during the August 3-4, 1999 
pBeppoSAX observation. Arrows indicate 500 s intervals for which only upper limits (3 a) can be set. The solid lines are the same as in Fig. 
hi The other curves show different applications of the stellar wind and centrifugal transition model presented in Section 4. The dashed, solid 
and dotted lines refer to angles x = 28°, 35°, 80°, respectively. These are the minimum angle compatible with the assumed mass inflow 
rate variation (28°), the maximum value (we take 80° since 90° provides the same result but more computational problems) and an angle 
providing a curve in between (35°). In the fit, ^ = 1 and the Be star equatorial disk model described in Section 4 have been assumed. 
Right panel: Pulsed fraction versus 0.1—200 keV unabsorbed luminosity during the BeppoSAX observations. Pulsed fractions were obtained 
by fitting the folded light curve with a sinusoidal signal. Error bars are a 1 cr. The curves represent the pulsed fraction k from the model 
presented in Section 4. No free parameter is involved except for the normalization fcmag (= 0.5). The dashed, solid and dotted lines refer to 
angles x = 28°, 35°, 80°, respectively. The upper limit (90%) refers to the BeppoSAX observation of Aug 13-16, 2000. 



accretion scenario (required by interpretation above), the 
issue as to why extreme variability manifests itself over 
a range of intermediate luminosities (between 2 x 10"^"' 
and 5 x 10'^^ ergs~^), while the source variability is much 
less pronounced for both higher and lower luminosities, 
remains unaddressed. 

The behaviour of 4U 0115+63 has a much more natural 
explanation in terms of the different regimes that are ex- 
pected for a rotating magnetic neutron star subject to a 
variable mass capture rate. The inflow of matter towards 
an accreting magnetic star is dominated by gravitational 
forces and possibly mediated by an accretion disk down 
to the magnetospheric boundary at r^. From this ra- 
dius matter can proceed inwards to the star surface only if 
the centrifugal force due to magnetospheric drag is weaker 
than gravity. This translates into the condition that the 
radius at which a test particle in Keplerian circular or- 
bit corotates with the central object [the so-called corota- 
tion radius, rcor = (G'MP^/4 7r^)^/^, with P the neutron 
star spin period and G the gravitational constant] is larger 
than Tni (lUarionov & Sunyaev 1975). In this direct accre- 
tion regime, the matter to radiation conversion efficiency is 
high L{R) = GM M/R and gives rise to luminous X-ray 
sources (here M is the mass inflow rate). If rm > rcor the 
drag by the rotating magnetosphere is super-Keplerian, 
such that the centrifugal force exceeds gravity and accre- 
tion is inhibited. However, as the inflowing matter reaches 
r„i, an accretion luminosity of L(rin) — GMAf/(2rin) 
must be released (the factor of 1/2 comes from assumption 
that the flow down to r^ is mediated by a disk). L{r^) 
is expected to scale approximately as '^ Af^/^. Therefore, 
above and below the transition to the propeller regime a 



smooth, close to linear, relationship between the accretion 
luminosity and M is expected. 

An important prediction is that the accretion luminosity 
across the transition from direct accretion to the propeller 
regime (corresponding to r-^ ~ ^cor) or vice versa should 
be characterized by a sudden luminosity jump 



iliin(-R)/^lim(rni) ^ 2 Tc 



oi/R 

/GMP^\ 



(1) 



1/3 



This is a factor of -^ 800 in the case of 4U 0115-^63 (Corbet 
1996; Campana & Stella 2000). In practice the transition 
separating the two regimes is expected to take place over 
a finite, though small, range of mass inflow rates around 
-^'^lim, such that a very steep dependence of the accretion 
luminosity on M ensues temporarily. 

Since the magnetospheric radius changes in response to 
variations of the mass inflow rate, the absolute luminosity 
at which the centrifugal barrier is expected to close can 
be determined based on models of the interaction between 
the inflowing matter and the rotating magnetosphere. In 
a simple spherical accretion approximation r^ is expected 
to scale as M~^/^ (Davidson & Ostriker 1973). A simi- 
lar dependence is obtained in the detailed models of the 
disk/magnetosphere that were developed by (e.g.) Ghosh 
& Lamb (1979) and Wang (1996) in the regime in which 
the disk is dominated by gas pressure (the one relevant 
to the case of 4U 0115-1-63). In order to account for the 
predictions of different models, we adopt the scaling above 
and introduce a correction factor ^, the ratio of the magne- 
tospheric radius determined on the basis of a given model 



to that of simple spherical accretion. We obtain 



Lii„,(i?) ^ 1 X 10 



37;r7/2 r2 






-7/3^^-2/3 



i?^ergs-i (2) 



where ^ = 1 for spherical accretion (by definition), ^ = 0.5 
in the model by Ghosh & Lamb (1979) and ^ ^ f in the 
model by Wang (1996). We take this range as representa- 
tive of the accuracy with which r-m can be predicted by 
current models (Mi. 4 and Rq are the neutron star mass 
and radius in units of 1.4 solar masses and 10^ cm, re- 
spectively. See e.g. Campana et al. 1998). Since the 
source distance (and, therefore, luminosity) and neutron 
star magnetic field are fairly accurately measured in the 
case of 4U 0115-^63, we obtain L\i^{R) ~ f^/^ 10^^ ergs'^ 
correspondingly Liim(rni) ~ ^^/^ 2 x 10'^'* ergs^^. The 
corresponding lines for ^ = 1 (continuous) and ^ = 0.5 
(dashed) are shown in Fig. n^. It is immediately appar- 
ent that the variations of the August 1999 observation fall 
just in the luminosity range of the transition between the 
propeller and the direct accretion regime, where very large 
luminosity variations are expected in response to modest 
changes of the mass inflow rate. On the contrary the qui- 
escent state luminosity of the August 2000 observation, 
during which the source flux was approximately constant, 
lies in the propeller regime, while the luminosity during the 
outburst observation on March 1999 is well in the range of 
the direct accretion regime. 

5. A SIMPLE MODEL FOR THE CENTRIFUGAL 
TRANSITION REGIME 

In this section we develop a simple model for the cen- 
trifugal transition regime and compare its predictions 
to the results from the August 1999 observation of 4U 
0115-1-63. In general, we express the source accretion lu- 
minosity in the transition regime as the sum of two con- 
tributions: (a) the luminosity of the disk extending down 
to the magnetospheric boundary, Ldisk! we assume that 
mass flows through the disk at rate M, which is equal to 
the rate at which mass is captured at the neutron star 
accretion radius; (b) the luminosity released within the 
magnetosphere imag by the fraction / of the mass inflow 
rate that effectively accretes onto the neutron star surface. 
We have 



L- 



^disk 



L„ 



L(r„0 + /(L(i?)-L(r„,)) (3) 



GMM [1/2 Tm + / (1/i? - 1/2 r„ 



This approximation the direct accretion regime corre- 
sponds to / = 1 and the propeller regime to / = 0. 

In the case of 4U 0115+63, being 2rcor -^ 800 i?, 
f L{R) > (1 - f)L{r^) for / > lO^^, i.e. the luminosity 
produced by matter accreting onto the neutron star sur- 
face is dominant over most of the centrifugal transition. 
This has two important implications. Firstly, being domi- 
nated by the release of energy within the magnetosphere, 
the emitted X-ray spectrum should be similar to that 
observed in the direct accretion regime (i.e. source out- 
bursts) for luminosities < 10'^^ ergs~^ and remain nearly 
unchanged across most of the centrifugal transition. We 
note also that according to models of accretion columns 
onto magnetic neutron stars, the spectrum is virtually in- 
sensitive to accretion rate variations as long as the optical 
depths remains < 1 (Nagel 1981). Only for luminosities 



of ^ lO"^"* ergs^^ (corresponding to / ^ 10^"^), when Ldisk 
becomes non- negligible, substantial spectral changes are to 
be expected. The apparent (relative) stability of the Bep- 
poSAX spectra during the August 1999 observation are 
consistent with this expectation. Secondly, as the transi- 
tion regime involves only a relatively small variation of ri„ 
the geometry of accretion close to the polar caps should 
show only minor changes (Wang & Welter 1981; Parmar 
et al. 1989). Therefore the pulse fraction is expected to 
remain essentially unaltered as long as Lmag dominates. 
On the other hand, Ldisk is expected to be unpulsed and 
the source pulsed fraction should decrease close to bottom 
of the transition regime. This is also consistent with the 
results from the August 1999 observation of 4U 0115+63 
(see Table 1 and Fig. |). 

We explored a simple model to determine / and the 
source behaviour across the centrifugal transition: we con- 
sidered a magnetic dipole field the axis of which is tilted 
relative to the neutron star rotation axis by an angle x- 
We determined the azimuthal ((/)-)dependence of the mag- 
netospheric boundary in the disk plane by equating the 
ram pressure of radially free-falling matter with the local 
magnetic pressure. The magnetospheric boundary takes 
an elongated shape given by 

r^{M,<t>)^^B^"R^^" (4) 

[1 + 3sin2 xsin^ 'l^?''^ (2 GA/M2)-i/7 

(Jetzer, Strassle & Straumann 1998). The minimum ra- 
dius rm(Af,0) corresponds also to x = 0; the approxima- 
tion usually adopted in models of the disk-magnetosphere 
interaction. The maximum radius rm(Aif, 7r/2) is only a 
factor of [1 + 3 sin^ xY'"^ < 1-49 larger. If rcor < ^^(M, 0) 
the magnetospheric boundary is larger than the corota- 
tion radius for any cj) and the propeller regime applies (i.e. 
/ = 0). If on the contrary Tcor > ''m(M,7r/2) every point 
on the boundary is within the corotation radius and the 
standard accretion regime applies (/ = 1). In the inter- 
mediate regime, rni(M,0) < rcor < ?'in(M,7r/2), only the 
fraction / ~ A(J)/2tt of the magnetospheric boundary for 
which Tcor > rin(M, 0) leads to accretion onto the neu- 
tron star surface and the transition regime applies. The 
accretion luminosity versus mass inflow rate curves have 
been calculated by using this model with ^ = 1. Different 
curves correspond to different values of the angle Xj the 
sharpest transitions occurring for the lowest values of x, 
as expected (see Fig. ||). 

We adopted the model above also to describe the evolu- 
tion of the source pulse fraction versus luminosity (see Fig. 
|2|). If Ldisk is unpulsed, the observed pulsed fraction k can 
be expressed as fc ~ fcmag f/{R/'2 ''m + /), where fcmag, the 
pulsed fraction of Lmag alone, can effectively be regarded 
as a free parameter, k depends on ^ (which sets the onset 
luminosity of the transition regime) and x (which deter- 
mines the range of the mass inflow rates over which the 
transition regime applies). Here too the dependence on x 
is such that for low values the transition is sharp, whereas 
for large values the transition occurs over a mass inflow 
rate variation by ~ 4 (see Fig. S). 

Modelling the time evolution of the source luminosity 
during the August 1999 observation would require know- 
ing the time evolution of M. We approximated M{t) with 
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Fig. 3. — Fraction of accreting mass / as a function of the mass inflow rate M normalized to the maximum mass inflow rate in the propeller 
regime, on the basis of the model discussed in Section 5. The different curves refer to angles x between 0° (left) to 90° (right) in step of 10°. 



the mass capture rate by the neutron star in its orbital 
motion within the equatorial disk of the Be star compan- 
ion over the phase interval 0.94-0.97. We assumed a 10 
km s~^ nearly Keplerian disk-outflow velocity and a den- 
sity profile oc r~'^ (as in Raguzova & Lipunov 1998). This 
model predicts a M variation of ;S 3 over the phase inter- 
val covered by the August 1999 observation. We remark 
that this is among the largest mass capture rate varia- 
tions that wind models predict. The centrifugal transition 
model described above with x — 35° and ^ = 1, was then 
used to calculate the line shown in Fig. ||; this reproduces 
reasonably well the overall behaviour of the light curve. 
Accordingly a luminosity variation of a factor of ^ 250 is 
produced in response to a factor of ~ 2 variation in M: 
clearly the centrifugal transition works like a very efficient 
amplifier. Moreover, these values of x ^^d ^ are among 
those that match nicely the pulsed fraction versus lumi- 
nosity variation (see Fig. ||) . By adopting the same model 
for the equatorial wind of the Be star over the phase inter- 
val of the August 2000 observation (0.42-0.51), the pre- 
dicted M is a factor of ~ 10 smaller (a value for which the 
neutron star is well in the propeller regime) giving rise to 
a luminosity in good agreement with the observed value. 

6. CONCLUSIONS 

The BeppoSAX observations of the transient X-ray pul- 
sar 4U 0115-1-63 revealed for the first time the presence of 
the extreme variations as the source approached perias- 



tron. The luminosities encompassed by these variations 
are within the range predicted by modelling of the cen- 
trifugal transition, which separates the propeller from the 
direct accretion regime over a small interval of mass inflow 
rates. On the contrary during a quiescent state observa- 
tion the source was most probably in the propeller regime, 
as its luminosity was lower still and nearly constant (how- 
ever see below). 

Before the present work the evidence for the onset of the 
centrifugal barrier was based on the sudden steepening of 
the outburst decay below a luminosity of ~ 10'^^ ergs"-^ in 
V0332-h53 (Stella et al. 1986) and 4U 0115-h63 (Tamura 
et al. 1992). In both cases, the source became quickly 
undetected in the relevant (collimator) instruments, such 
that the transition towards the propeller regime could not 
be seen. It was also proposed that several other sources 
in their quiescent (or low) state host a neutron star in the 
propeller regime: among these are the HXRTs A 0538- 
66 (Campana 1997; Corbet et al. 1997) and A 0535-^26 
(Negueruela et al. 2000). The evidence reported here is 
far more convincing in that dramatic source flux variations 
of the kind expected in the centrifugal transition regime 
were observed for the first time. These variations occur 
within the luminosity interval predicted by models of disk- 
magnetospheric interaction for the case of 4U 0115-1-63, 
where the neutron star spin and magnetic field, as well as 
distance are fairly accurately measured. The source flux 
level and absence of sizeable variability in the faint state 



away from periastron is also in agreement with basic ex- 
pectations for the regime in which the centrifugal barrier 
is fully closed. Yet, owing to poor statistics a thermal 
like spectrum cannot be ruled out. The analogy with neu- 
tron star soft X-ray transients suggests that such a spec- 
tral component, if present, might be due to reemission of 
heat in the inner crust caused by pycnonuclear reactions 
(Brown, Bilsten & Rutledge 1998; Campana et al. 1998; 
Colpi et al. 2001). The inferred blackbody luminosity^ is 
in ~ lO'^^ ergs~^ to be compared with an expected deep 
crustal heating luminosity of ~ 5 x 10'^'^ ergs~^ under the 
hypothesis that 4U 0115-1-63 has accreted at a time average 



rate of ^ 5 x lO^'^ gs~^ (as suggested by the RossiX TE 
ASM light curves, b.ttp://xte.mit .edu/lcextrct/ ) for 
some ~ 10^ years. Independent of the origin of the qui- 
escent emission from 4U 0115+63, the data presented in 
this paper provide substantial new evidence in favor of the 
centrifugal transition regime. 
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